Excitation Spectrum at the Yang- Lee Edge 
Singularity of the 2D Ising model 



TOMASZ WYDRO 

Laboratoire de Physique Moleculaire et des Collisions, Universite de Metz, 
1 bvd Arago, 57078 Metz, France 
wydro@sciences. univ-metz . fr 

JOHN F. McCABE 
412 Morris Ave., # 34, Summit, NJ 07901, USA 
jfmccabe2@earthlink.net 

February 2, 2008 



Abstract 

This paper studies the Yang-Lee edge singularity of 2-dimensional 
(2D) Ising model through a quantum spin chain. In particular, finite-size 
scaling measurements on the quantum spin chain are used to determine 
the low-lying excitation spectrum and central charge at the Yang-Lee edge 
singularity. The measured values are consistent with predictions for the 
(AijAi) minimal conformal field theory. 



In 1978, Fisher[l] argued forcefully that Yang-Lee edge singularities [2, 3] 
are similar to ordinary critical points. Later, Cardy[4] provided another strong 
argument that the Yang-Lee edge singularity of the 2D Ising model is identi- 
fied with the (A 4 ,A 1 ) minimal conformal field theory (CFT)[5, 6] of the ADE 
classification. [7] Cardy's identification enables a number of CFT predictions 
about the Yang-Lee edge singularity of the 2D Ising model. 

Some of these CFT predictions have been confirmed. Using transfer matrix 
methods, Itzykson et al measured the central charge and exponent v at the 
Yang-Lee edge singularity of the 2D Ising model. [8] Using short quantum spin 
chains, [9] Uzelac et al also measured the exponent v at the Yang-Lee edge 
singularity of the 2D Ising model. [10] Both measurements arc in agreement 
with predictions derived from Cardy's identification. 

This article provides other measurements that can be predicted based on 
Cardy's identification. The measurements provide the low-lying excitation spec- 
trum at this Yang-Lee edge singularity of the 2D Ising model. This article also 
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compares the measured low-lying excitation spectrum with predictions based 
on Cardy's identification of the {A±,A{) minimal CFT with this Yang-Lee edge 
singularity of the 2D Ising model. [4, 8, 10] 

In this article, measurements]!!] are based on finite-size scaling in quantum 
spin chains. For the 2D Ising model in an imaginary external magnetic field, the 
associated iV-site quantum spin chain has a Hamiltonian, Hi S i ng , given by: [12] 

N 

Hi S i n g = - '^2{ta z (n)<r z (n + 1) + iBcr z (n) + a x (n)} . (1) 
n=l 

In Eq. (1), a x (n) and a z (n) are 2x2 Pauli spin matrices at site n, the parameter 
"i" is a positive coupling for ferromagnetic spin-spin interactions, and iB is a 
purely imaginary external magnetic field. In Eq. (1), the last term results from 
single inter- row spin flips in the associated 2D transfer matrix. [13] 

In this article, the phcnomcnological renormalization group (PRG) is used to 
determine critical values of the purely imaginary magnetic field, iByl{N) 1 for 
various chain lengths, N. For such purely imaginary magnetic fields, the PRG 
equation becomes: [14, 8] 

[N - l]m(B YL {N),N - 1) = [N]m(B YL (N),N) . (2) 

In Eq. (2), m(B,N) = [E^B,^ -E (B,N)] where E (B,N) and E\{B, N) 
are the energies of the ground state "0" and first excited state "1" of the quan- 
tum spin chain having length N. At the Byl(NYs, the Ising quantum spin chain 
exhibits finite-size scaling behavior associated with the Yang-Lee edge singular- 
ity of the 2D Ising spin model. In particular, if the Byl(N)'s converge to a 
nonzero value as N — > oo, that value corresponds to the critical point for the 
Yang-Lee edge singularity of the 2D Ising spin model. 

At the Byl{NYs, excitation energies and other physical quantities should 
scale. CFT predicts the scaling behavior of such physical quantities at the 
B Y l(NYs. 

First, CFT predicts how excitation energies will scale with the length, N, of 
the quantum spin chain. For an excited energy eigenstate "i" of the quantum 
spin chain, the CFT prediction is that the excitation energy, Ei(N) — Eo(N), 
will scale as: [15] 

E l (N ) -E (N)^C2 ^ + ^ N {A + A) ■ (3) 

In Eq. (3), Aj and Aj are left and right conformal dimensions of the field 
"i" in the associated CFT, and A and A are the conformal dimensions of the 
primary field of lowest "negative" scaling dimension in the CFT. Such negative 
dimension fields occur in various non- unitary CFTs. In Eq. (3), the constant ( 
is non-universal and depends, e.g., on the normalization of the Hamiltonian of 
the quantum spin chain. 

Second, CFT predicts how the ground state energy, Eq(N), will scale with 
the length, N, of the quantum spin chain. In particular, the ground state energy 
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scales as: [17, 16] 

E (N) = AN + B-(irC eff /(6N) + ... (4) 

While the constants A and B are non-universal, C e // and the exponents of 
the finite-size corrections in 1/N of the last two terms are universal numbers 
predicted by CFT. In particular, C e ff, is the effective central charge of the 
CFT. Ceff is defined by C e ff = C - (A + A) where "C" is the normal central 
charge of the CFT. [8] 

For the minimal CFTs, the ADE classification provides modular invariant 
partition functions [7, 1 8] from which the low- lying excitation spectrum and cen- 
tral charge are easily extracted. For the (Ai,A\) minimal CFT, Table 1 shows 
the low-lying excitation energies and associated degeneracies of the spectrum 
obtained from the associated modular invariant partition function of the ADE 
classification. Rather than absolute excitation energies, Table 1 lists normal- 



CFT 


{A i ,A 1 ) 




Normalized Energies 





1 


2.5 


5.0 


6.0 


7.5 


Degeneracy 


1 


1 


2 


3 


2 


4 



Table 1: Low- lying excitation spectrum of {A±,A{) minimal CFT 

ized excitation energies, i.e., ratios. For a state "i", the associated normalized 
excitation energy is defined as the ratio is the actual excitation energy of the 
state "i" over the actual excitation energy of the lowest excited state " 1" . Here, 
absolute excited energies are measured with respect to the ground state energy. 
The ratios of Table 1 have the advantage of not depending on non-universal 
constants such as £. For the {A±,A\) minimal CFT, the form of the modular 
invariant partition function also implies that the effective central charge, C e t / , 
is equal to 2/5. [8] Below, finite-size scaling measurements on the Ising quantum 
spin chain are used to determine the low-lying excitation spectrum and effective 
central charge at the Yang-Lee edge singularity. 

The measurements of the critical magnetic fields, Byl(N), were obtained by 
numerically solving PRC eq. (2) for Ising quantum chains of different lengths. 
For the numerical solutions, the state energies were obtained by applying the 
Lanczos algorithm to Hismg of eq. (1). Table 2 shows the values for the critical 
fields, i.e., Byl(NYs, ground state energies, and lowest excitation energies, i.e., 
Gap(N)'s, as obtained via the Lanczos algorithm. Table 2 lists measurements of 
these physical quantities for Ising quantum spin chains in which the coupling, 
t, had the value of 0.1. The measurements of Table 2 also show that NxGap(N) 
scales to a constant as N — > oo as expected for the PRC 

Based on the finite-size scaling measurements of Table 2, we first evaluated 
the effective central charge, C e //, at the Yang-Lee edge singularity of the 2D 
Ising model to test the numerical methods. In particular, ground state energies 
for adjacent triplets of chain lengths [(N-l), N, (N+l)] and an estimate, ((N), 
for the non-universal constant, (, were used to solve Eq. (4) for estimates, i.e., 
C e ff(N)'s, to the effective central charge, C e //. In each evaluation, the Gap(N) 
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l\Tn m nPT" 

of Sites 


Byl(N) 


Energy of ground state 


Gap(N) 


N x Gap(N) 


3 


.2459180i 


-2.8811043 


.8103423 


2.4310 


4 


.2384127i 


-3.8028211 


.6629112 


2.6516 


5 


.2352339i 


-4.7341982 


.5613016 


2.8065 


6 


.2337637i 


-5.6688215 


.4858628 


2.9152 


7 


.2330279i 


-6.6048003 


.4275400 


2.9928 


8 


.2326347i 


-7.5414746 


.3811698 


3.0494 


9 


.2324118i 


-8.4785910 


.3435105 


3.0916 


10 


.2322793i 


-9.4160213 


.3123765 


3.1237 


11 


.2321972i 


-10.353696 


.286250 


3.1488 


12 


.2321442i 


-11.291568 


.264041 


3.1685 


00 


.23193i 


— oo 


0.0 


3.2840 



Table 2: PRG results for Byl{N) 1 Ground state energy, Gap, and NxGap as a 
function of the number of sites, N 



was used to estimate the non-universal constant ((N) of Eq. (3) at chain length 
"N". The resulting estimates, i.e., the C e //(./V)'s, are plotted in Figure 1. The 
C e f f (N) 's have a " 1 /N" dependence due to the higher order " 1 /N" corrections 
to eqs. (3) and (4). 

From a visual inspection of Figure 1, one sees that the estimates for the 
effective central charge have a crossover behavior in 1/N. After the crossover, 
the C e //(iV)'s converge smoothly towards a limit, i.e., C e //, as 1/N — > 0. The 
effective central charge may be obtained by extrapolating the estimates towards 
1/N = 0. A visual inspection of Figure 1 also shows that the C e ff(N)'s behave 
as C e ff + a/N as 1/N — ► 0. Based on this form, one easily extracts from 
the measured values that C e ff is approximately equal to 0.41. This measured 
value for C e // agrees well with the value of 2/5 that the (At, A\) minimal CFT 
predicts for C e //. 

The critical magnetic field values of Table 2, i.e., the Byl's, were also used 
to find the low-lying excitation spectra of Ising quantum spin chains of vari- 
ous lengths. While the Lanczos algorithm can provide the low-lying excitation 
spectra, it is inconvenient for determining the complete low-lying excitation 
spectrum. In particular, the Lanczos algorithm requires one to select a starting 
vector and then, produce a Krylov space over the starting vector. For simple 
choices of the starting vector, we found that the Lanczos algorithm spreads 
the low-lying excitations over several Krylov spaces. Thus, an evaluation of the 
complete low-lying excitation spectrum via the Lanczos algorithm would require 
finding a first Krylov space and then, finding other Krylov spaces orthogonal to 
the first space. To avoid such numerical complications, we used the commercial 
linear algebra package of MaplerM 9.5 to generate complete excitation spectra. 

The MapleTM package provided excitation spectra for Ising quantum spin 
chains with 6-12 sites. The measured low-lying parts of said spectra including 
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Figure 1: The C e ff(N) estimates for the effective central charge, C e //, as a 
function of 1/N. 



both energies and degeneracies are provided in Table 3. There, the excitation 
energies are again normalized by division by the lowest excitation energy of the 
Ising quantum spin chain. As already described, such a normalization removes 
any dependence on the non-universal constant 

From the measured spectra of Table 3, one can determine the form of the 
low-lying excitation spectrum in the limit where 1/N —> 0. The determination 
requires first, identifying classes of eigenstates that correspond for different chain 
lengths and then, determining the scaling behavior for each of the identified 
classes as, 1/N — > 0. The identification of corresponding eigenstates for different 
values of N is achieved with the aid of some simple rules. The first rule is that 
the energies of corresponding states vary monotonically and smoothly with N 
provided that the initial value of N is sufficiently large. The second rule is 
that the correspondences between states must account for degeneracies. In 
particular, while new eigenstates appear as N increases, old eigenstates do not 
disappear. 1 Thus, in each class, degeneracies of eigenstates either stay constant 
or increase with N. 

1 A few "special" low- lying excited eigenstates have excitation energies whose magnitudes 
grow rapidly with N. The" special" eigenstates are not part of the low-energy excitation spec- 
trum as N — > oo. In fact, the number eigenstates with energies of smaller magnitude than 
those of the "special" eigenstates grows with N. Such "special" states were also seen at the 
Yang-Lee edge singularity of the 3-state Potts quantum spin chain where they were not part 
of the spectrum in the thermodynamic limit. [21] The "special" excited eigenstates are not 
discussed further. 
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State /[Degeneracy] 


6 


7 


8 


9 


10 


11 


12 


A/ [2] 


2.68432 


2.64386 


2.61415 


2.59207 


2.57540 


2.56260 


2.55253 


B/[l] 


4.18193 


4.27896 


4.36713 


4.44474 


4.51197 


4.56977 


4.61912 


C/[2] 


4.51738 


4.63236 


4.70368 


4.75182 


4.78652 


4.81281 


4.83329 


D/[2] 


5.85889 


5.89208 


5.91240 


5.92644 


5.93703 


5.94544 


5.95210 


E/[2] 




5.68559 


6.03104 


6.27270 


6.45018 


6.58573 


6.69223 


F/[2] 




6.24252 


6.35798 


6.46344 


6.55966 


6.64694 


6.72535 



Table 3: The table shows normalized low-lying excitation energies of states A - 
F for Ising quantum spin chains with 6 to 12 sites. Degeneracies of the states 
are as listed in the left column. 



Application of the above rules leads the classification of low-lying excited 
eigenstates of Table 3. The classification includes classes A - F of state types. 2 
Figures 2-5 plot the normalized excitation energies for the states of classes A 
- F as a function of the inverse of the length of the Ising quantum spin chain. 
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Figure 2: The measured energies of type A states are plotted in 1/N. 



A visual inspection of Figures 2-5 readily shows that the excited states A, 
B, C, D, E, and F combine into the four distinct sets A, B & C, D, and E & 
F. Within each distinct set, the excited eigenstates have energies that approach 
the same value as 1/N — > 0. 

2 Here, the lowest excited state has been ignored. The lowest excited eigenstate was how- 
ever, observed to be non-degenerate as predicted for the (A^,A\) minimal CFT. The lowest 
lying excited eigenstate automatically has an energy of 1 in the convention of Table 1. 
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Figure 3: The measured energies of type B states (squares) and of type C states 
(circles) are plotted in 1/N. 
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Figure 4: The measured energies of type D states arc plotted in 1/N. 
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Figure 5: The measured energies of type E states (squares) and of type F states 
(circles) are plotted in 1/N. 



A visual inspection of Figures 2-5 shows that the normalized energies of 
the states of the sets A, B & C, D, and E & F approach about 2.45, 5.0, 6.03, 
and 7.6, respectively, as 1/N — > 0. A BST analysis[20] was performed to obtain 
more reliable values for the energies of these states in the limit where 1/N — > 0. 
The BST analysis indicated that the normalized energies of the states of type 
A, B, C, D, E, and F scale to 2.4995(5), 5.005(1), 5.003(3), 5.99(1), 7.54(8), 
and 7.60(7), respectively, in this limit. 3 These finite-size scaling measurements 
of the low- lying excitation energies agree well with those of the (A4, A) CFT 
as shown in Table 1. 

Also, an inspection of Table 3 shows that the distinct state sets A, B & 
C, D, and E & F, of different limiting excitation energies have 2, 3, 2, and 4 
states, respectively. Thus, the PRG measurements of the normalized excitation 
spectra also provide values for the degeneracies that agree well with those of 
the {A^A-l) CFT of Table 1. 

In conclusion, our finite-size scaling measurements on the Ising quantum 
spin chain at the Yang-Lee edge singularity have produced a low-lying excitation 
spectrum that is in very good agreement with that predicted from the (Aj, A{) 
minimal CFT. These results further confirm Cardy's identification of the Yang- 
Lee edge singularity of the 2D Ising model with the (At, A\) minimal CFT. 



3 The BST analysis also provided evidence that finite-size scaling corrections to the nor- 
malized energies of these sets of states scale as 1/N b where the b's are in [1, 2] as 1/N — > 0. 
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